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I define a statistical model of graphs in which 2-dimensional spaces arise at low 
temperature. The configurations are given by graphs with a fixed number of edges 
and the Hamiltonian is a simple, local function of the graphs. Simulations show that 
there is a transition between a low-temperature regime in which the graphs form 
triangulations of 2-dimensional surfaces and a high-temperature regime, where the 
surfaces disappear. I use data for the specific heat and other observables to discuss 
whether this is a phase transition. The surface states are analyzed with regard to 
. topology and defects. 
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I. INTRODUCTION 

O'. 

Spacetime represents one of the remaining frontiers in theoretical physics. We have 
successful theories of fields in spacetime and of the geometry of spacetime. But there is no 
established theory of spacetime itself, a theory that would explain why we live in a spacetime 
and not in something completely different. 
■ There are approaches that go in this direction, notably string theory, matrix models 

0-0, twistor theory |5|, ||, causal dynamical triangulations |7|, [|, Regge calculus ||, loop 
quantum gravity [0-[12[ and quantum graphity Some of these start from objects 

that are radically different from spacetime and one hopes that spacetime arises from the 
dynamics of these objects. In others, spacetime is assumed from the beginning and only the 
topology and geometry is dynamical. 

In this article, I am concerned with the former type of approach, i.e. the attempt to obtain 
space from a theory, where there is no space to start with. I define a simple statistical model 
of graphs that is designed to have 2-dimensional triangulations as ground states. The setup 
is close in spirit to earlier models of quantum graphity explored by Konopka, Markopoulou, 
Smolin and others [[H|-[[rj|. I present data from Monte Carlo simulations 1 of graphs with 
45, 90 and 180 edges. The results show that the system settles into 2d triangulations at 
low temperatures. Above a transition temperature the surface states disappear and more 
general graphs dominate. As far as I know, this is the first time that a model of graphs 
succeeds in reaching manifold-like states without help from additional ad hoc constraints. 
The temperature dependence of specific heat and other observables exhibits similar features 
as in phase transitions of more conventional systems. The transition temperature, however, 
drifts as the size of the system is changed, leaving both a finite- and zero-temperature phase 
transition as a possible scenario for the infinite-size limit. 

In a second set of simulations, I investigated the surface properties in more detail with 
regard to topology, defects and curvature. At very low temperature the surfaces are typically 



5h . 



X 

5-( 



Electronic address: fconrady@perimeterinstitute.ca 



1 For an introduction to Monte Carlo methods, see e.g. p 
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connected and non-orient able, and carry a number of conical singularities. 

In the space of graphs the surface states form a very complex subset. Like in other 
systems with a complex energy landscape (see e.g. [p0|-p2|), this makes it difficult to ensure 
equilibrium at low temperature. The simulation may get trapped in energy valleys and 
sample no longer according to the Boltzmann distribution. To check whether the simulation 
is equilibrated, I compare the original simulation with an alternate simulation schedule that 
involves frequent restarts and annealing. 

Although triangulations play a role in this model, the setup is quite different from causal 
dynamical triangualations and Regge calculus. In these theories, each individual configu- 
ration is a space in the form of a triangulation, and one finds different regimes, where the 
effective geometry is smooth, crumbled or polymer-like, for example. In our case, the con- 
figurations are graphs and the regimes are distinguished by the existence or absence of a 
space. 

For testing and inspecting the simulation I used the Ubigraph software for visualizing 
dynamic graphs p3 |. Videos of the simulation can be viewed at this link f2~4|| . 

The paper is organized as follows. In section I define the model and discuss the 
expected ground states and its defects. I also define an observable used to measure the 
"superficiality" of graphs. Section |I| describes the details of the algorithm — the Monte 
Carlo moves and the trial probabilities needed for detailed balance. This is followed by sec. 
[y] , where I report the results of the simulation. I discuss evidence for and against a phase 
transition and the graphs' properties at low temperature. The question of low temperature 
equilibration is scrutinized in sec. |V|. In the final section, I give a summary and discuss some 
of the conceptual aspects of the model. 



II. THE MODEL 
A. Definition 

The model is a statistical model of graphs. Configurations are given by undirected, 
labelled graphs T with a fixed number N e of edges. The edges of the graphs are labelled by 
numbers 1,2,. .. , N e . Automorphic graphs are distinguished if they differ in their labelling 
of edges. Vertices have to be contained in at least one edge. There can be at most one edge 
between two vertices, i.e. double, triple etc. edges are not allowed. Moreover, edges from a 
vertex to itself are prohibited. 

Vertices and edges of graphs are denoted by v and e respectively. When n + 1 vertices 
form a complete subgraph (i.e. when every vertex of the subgraph is connected with each 
of its other vertices), I call this subgraph an n-simplex. Thus, vertices and edges define 
0- and 1-simplices. 2- and 3-simplices will be referred to as triangles t and tetrahedra r 
respectively. For the definition of the Hamiltonian the notion of valence is important. An 
n-simplex A is said to have valence Va, if it is contained in Va (n + l)-simplices. 

The Hamiltonian is specified in such a way that it favors graphs corresponding to a 2- 
dimensional triangulation. That is, graphs are energetically preferred when its subgraphs 
form simplices of a 2-dimensional simplicial complex. The Hamiltonian depends on the 
graph T and is a sum of four types of terms that are associated to vertices, edges, triangles 
and tetrahedra: 

v e t t 
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The sums extend over all vertices v, edges e, triangles t and tetrahedra r contained in T. 
The individual terms are defined as follows: 

( , V e < 2 , 

H v = c v \V V - 6| , H e = < , H t = -c t} H T = c T . (2) 

[ Cf y e 2 , k>2, 

The coefficients c v , c e , c t and c r are positive constants and their specific values are taken to 
be 

c v = 1 , c e = 10 , c t = 7 , c T = 500 . (3) 

The most important terms are H e and H t . The triangle term H t favors the creation of 
triangles, since every triangle lowers the energy by c t . The edge term H e penalizes edges with 
a valence V e greater than 2. It thus suppresses the branching of triangles, which would be 
incompatible with a 2d triangulation. The tetrahedral term H T suppresses the formation of 
tetrahedra and allows us to exclude surfaces that consist merely of disconnected tetrahedra. 
The vertex term H v was introduced in order to favor smooth, low-curvature triangulations 
over rough, and crumbled ones 2 . The curvature at a vertex v is determined by its valence 
V v , and in 2 dimensions a valence of 6 corresponds to flatness. 
The partition function is defined by 

Z = J2e~^, (4) 
r 

where /3 is the inverse of the temperature T. The sum includes all graphs with N e edges 
that are admissible according to the rules stated previously. The average of an observable 
Or is given by 

<0>=4l> e " /Mrr - (5) 
r 



B. Surfaces, defects and their observables 

The above choice of the Hamiltonian makes it likely that triangulations of surfaces are 
among the ground states of the system. Surfaces, in the strict sense, need not be the only 
minima of energy, however, since there can occur defects to which the Hamiltonian is not 
sensitive. 

Through the identification of complete subgraphs with simplices, each admissible graph 
T defines a simplicial complex. A simplicial complex is a triangulation of a surface (with 
boundary) if and only if 3 

(a) all edges have valence V e = 1 or V e = 2, and 

(b) at every vertex edges and triangles form a sequence ei t\ &3 ■ ■ ■ e n t n e n+ i, where 
ej, e i+ i C U (and e x = e n+x or e x ^ e n+1 ). 



2 As such a triangulation defines only a topological manifold, but if we think of each edge as having the 
length 1, it becomes a manifold with a metric, i.e. a manifold with a notion of distance and curvature. 

3 see e.g. |§, sec. 1.4.22 
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The term H e in the Hamiltonian ([I]) does not distinguish between valence V e = 1,2 and 
V e = 0, so there could be minima that consist of surfaces dressed with additional edges 
outside of any triangle (see Fig. [I]). The term H t tends to eliminate these edges in order to 
maximize the number of triangles for a given total edge number N e . Nevertheless there could 
remain edges as "odd men out", since their inclusion would require a restructuring of the 
triangulation. The analysis is further complicated by the term H v , which is sensitive to the 
valence of vertices. Finally, even when condition (a) is fulfilled, the simplicial complex may 
not be a surface, because condition (b) is violated at one or several vertices. Such vertices 
correspond to so-called conical singularities, where two surfaces are joined in a single point 
(see Fig. 0). 

For these reasons, it is difficult to state precisely which graphs are the minima of energy. 
The previous arguments suggest that the ground state is given by surfaces with defects in 
the form of conical singularities and edges with valence V e = 0. 

The purpose of the present work is to determine by simulation whether this intuitive 
picture is correct and if the system does indeed settle into surface-like states at low temper- 
atures. Moreover, this should allow one to see how the system departs from these states as 
the temperature is increased. To quantify the presence (or absence) of a surface, I use an 
observable M that measures how close the graph is to being a triangulation of a 2d manifold: 
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o , v e = 



1/2, V e = l, 

(6) 



1 , V e = 2 



, V e > 2 . 



Edges with valence V e — 2 contribute with 1 to M, boundary edges add 1/2 and any other 
edge counts as zero. When comparing systems of different size N e , it is convenient to use 
the "intensive" quantity 

m = —M. (7) 

A graph with m = 1 is, up to conical singularities, a perfect surface without boundary. As 
m is not sensitive to conical singularities, the number iV con of such singularities is another 
observable that will be measured in the simulation. By analogy to the susceptibility of a 
magnet, I will also define and measure a "susceptibility" related to M: 

X = ^- {(M 2 ) - (M) 2 ) = /3iV e ((m 2 ) - (m) 2 ) (8) 
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III. MONTE CARLO SIMULATION 

I simulated the model by the Monte Carlo method using a Metropolis algorithm. In 
setting up this simulation the main challenge was to find suitable moves that allow one 
to explore the states of the system efficiently at low temperatures. One can adopt a very 
simple set of moves that displace edges randomly from one location to another, and such 
moves would be clearly ergodic. It is far from clear, however, that these moves are also 
powerful enough to sample all the relevant configurations within a limited simulation time. 
The surface configurations deemed to be important at low temperature constitute only a 
tiny fraction of the full space of graphs. It is doubtful that a mere hit and miss method 
can reach these very special graphs. Initial simulations showed, in fact, that such minimal 
moves are not sufficient for finding these configurations. 

To resolve this problem I extended this set of moves by more elaborate moves that 
assist the system in building low-energy configurations. The resulting algorithm succeeds 
in sampling the surface graphs 4 . 

In the next two subsections, I describe the moves of this algorithm and also the trial 
probabilities that are needed to achieve detailed balance. 

A. Moves 

Since the number N e of edges is conserved, each Monte Carlo move consists of a succession 
of two more basic moves. The first move removes m edges from the graph and the second 
one adds again m edges. 

I found it convenient to categorize these moves into 6 types, labelled by a pair of numbers 
n m . The number m denotes the change in the number of edges, while n labels the specific 
way in which the addition (or removal) of edges happens (see Fig. |2|). 

Roughly speaking, n m , m > 0, stands for a move that creates a new simplex by adding 
m edges. If n = 1, this is done by adding a new simplex that is disconncted from the rest 
of the graph. If n — 2, m edges are "glued" onto an existing (m — l)-simplex A m _! to form 
an m-simplex A m . If n = 3, the m-simplex A m is built by "glueing" m edges to an existing 
(m — l)-simplex A m _i and a vertex. 

More concretely, this means that the move li adds an edge whose vertices have both 
valence 1. The move 2i adds an edge such that one vertex of the edge has valence 1 and the 
other vertex has valence greater than 1. The 3i move adds an edge between two existing 
vertices. 

I also introduce a fourth move, called 4 X . It adds an edge between two vertices that are 
linked to a common vertex, but are not linked directly. These moves are redundant in the 
sense that they are equivalent to a subset of the 3i moves. The 4j moves are nevertheless 
important, since the explicit proposal of such moves makes the algorithm far more efficient in 
finding minimum energy configuration. They appear to facilitate the growth of 2-simplices 
and surfaces. 

In the case m = 2, we have, in analogy to 2i, the 22 move which "glues" two edges a and 
b to an existing edge c such that a, b, c form a 2-simplex. The vertex shared by a and b has 



4 It is another question whether this sampling is also unbiased, i.e. representative of the low-temperature 
equilibrium. This will be discussed in section M. 
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Figure 2: List of moves from which Monte Carlo moves are built. 

valence 2. The move 32 does the same as 2 2 except that the vertex common to a and b has 
valence greater than 2 (i.e. one end of a and b is "glued" to an existing vertex of the graph). 
Like the A 1 move, the 2 1 and 2 2 move are not needed for ergodicity, but they enhance the 
algorithm's efficiency at low temperature. 

In principle, one could also define a 1 2 move that adds a disconnected 2-simplex. However, 
I found it easier to omit this move in the algorithm. Since the edge number is conserved, 
1 2 would have to be accompanied by two moves of m — —1, which would require additional 
program code. 

The moves n m with negative m are defined as the inverses of those with m > 0. Altogether 
we have the following set of basic moves 5 : 

{l±i, 2±i, 3±i, 4x, 2 ± 2, 3 ±2 } 

Each of the Monte Carlo moves is given by a sequence of moves (n- m ,n' m ), m > 0, where 
n_ m and n' m are taken from the previous set. 

B. Trial probabilities and detailed balance 

Next I describe the way Monte Carlo moves are proposed. The resulting proposal (or trial) 
probabilities are generally not symmetric w.r.t. inversion of the move. Therefore, to ensure 
detailed balance, we need to include trial probabilities when determining the acceptance 
probability. 

The trial probability for the Monte Carlo move (n- m: n' m ) is given by 

P t (n- m ,riJ = P t (n- m )P t (ri m ), (9) 

where P 4 (n_ m ) and P t (n' m ) are the trial probabilities for the subtraction and subsequent 
addition of edges. 



5 We do not need 4_i, since 2-simplices are easily destructed by 3_i. 
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For the subtraction move n_ m , m > 0, the selection proceeds as follows. There is a 
probability P m of choosing m = 1 or m = 2 that is weighted in a specific manner by 
the number of existing simplices. Then, we have a probability P Am for choosing the m- 
simplex A m on which the move operates. This simplex is chosen randomly from the set of 
m-simplices, so 

p ^ = it- (10) 

Finally, one has to select the type n of the move. The probability P n for this is 1 in all cases 
except for 3_2 ; where it is P3 = 1/3, since there are three ways of removing two edges from 
a given 2-simplex. The overall trial probability for the subtraction move n_ m is 

P t (n^ m ) = P m P Am P n . (11) 

In the case of the additive move, we do not have to choose m, since it is fixed by the previous 
removal of edges. The choice of type n! is based on a certain weighting P n > whose details we 
omit. The non-trivial part is the probability P Tec (n' m ) for selecting the simplices that serve 
as "receptors" for the edges to be added. For moves of type 1, P rec (l m ) = 1. For n' = 2, we 
have ^ 

P rec (2 m ) = PA m _i = 77 , (12) 



m — 1 



where P\ m _i is the likelihood of picking one of the (m — l)-simplices at random. For n' = 3, 
one has to choose both an (m — l)-simplex and a vertex, therefore 

Prec(3 m ) = PAm-l^V = TT ■ (13) 

Am— l v 

The total trial probability for n' m is 

P t (n'J = P n ,P Tec (n'J . (14) 

An exception to this rule occurs when we deal with a move that is both 3i and 4\. In this 
case, we have 

P t (3i) = P t (40 = P 3 P rcc (3 1 ) + P 4 P rec (4i) , (15) 

where P rec (3i) is given as before and P rec (4!) is the likelihood of proposing the receptors of 
the 4i move. The latter are determined by first choosing a vertex v c that I call the "center" 
vertex and then vertices v\ and v-i that are linked to v c , but not linked to each other (see 
Fig. 0). Thus, the likelihood for a given pair (t>i,t>2) can be written as 

f -w4£r (16) 

Here, the sum extends over all center vertices v c compatible with the pair (v%, V2), while N Vc 
is the number of all pairs (1)1,1)2) for which v c is a center vertex. 

Once the trial probabilities are determined, we can evaluate the acceptance probability 
P a (n„ m ,n' m ) via 

P„( n _ m ,<) = mm {l , g("'-.." r ) -( g - E )/Tj . (17) 



E and E' denote the energy before and after the move. 
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Table I: Equilibration times and total number of Monte Carlo steps. 



IV. RESULTS OF THE SIMULATION 



In this section, I present the results of the simulation. Measurements were done for three 
sizes of the system, N e = 45, 90 and 180, and for a range of temperatures from T = 0.5 to 
T = 6.75. I begin with the estimates for equilibration times (the next subsection). Then, in 
subsection |IV B| , I describe the results for the two main observables, the energy E and the 
observable M (defined in eq. (§)), and the quantities derived from them, the specific heat 
per edge c and the susceptibility x associated to M (see eq. @). In subsection |IV Q , the 
low-temperature regime is investigated in more detail with regard to topology, curvature 
and defects. In section |V|, I will discuss the issue of equilibrium at low temperature and 
caveats entailed by it. 



A. Equilibration times 

Figure |] contains examples of plots used to estimate the equilibration times for different 
temperatures T and sizes N e of the system. The plots show data from preparatory simu- 
lations that measured the energy per edge e = E/N e and the observable m = M/N e for 
a duration of t — 2 • 10 7 Monte Carlo steps. The starting configuration is a graph of N e 
disconnected edges 

It should be noted how the characteristics of the plots change as we go from low to higher 
temperatures. The plot is very autocorrelated for T = 0.5, becomes less correlated and more 
fluctuating at T = 1.125, and turns into strong fluctuations with little autocorrelation at 
T = 1.5. 

Based on these plots and similar ones for other temperatures and edge numbers, I adopted 
a set of equilibration times At eq for the main simulation which are listed in table |. The 
total number of Monte Carlo steps t MC is taken to be 100At eq in each case. 
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Figure 3: Examples of plots that were used to estimate equilibration times t eq : energy per edge e 
and m versus simulation time t for N e = 180 edges. 

B. Evidence for a phase transition? 

Figure f| displays the simulation results for the energy per edge e and the observable m. 
Here and in the following error bars are omitted whenever they are smaller than the symbols 
of the data points 6 At sufficiently low temperature m reaches the value 1 and, up to conical 



Correlation times and error estimates will be addressed in section 
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Figure 4: Energy per edge e and m versus temperature T for edge number N e = 45, 90 and 180. 



singularities, the average graph is given by a surface without boundary. This confirms the 
expectations expressed in the previous section. The attainment of m = 1 goes along with the 
energy per edge reaching its bottom value. (Curiously, for N e = 45, the energy rises again 
slightly as T decreases further. This might be related to incomplete equilibration at low 
temperatures, see subsec. 0). The transition temperature, where m — 1 sets in, decreases 
as N e increases. 

That something interesting happens when m becomes 1 is also suggested by the results 
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Figure 5: Specific heat per edge c versus temperature T for N e = 45, 90 and 180. 
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Figure 6: Specific heat per edge c versus temperature T. 



for the specific heat, shown in Fig. ||. The specific heat can, in principle, be determined in 
two ways; from the variance, 

c = f? Ne « e 2) - ( e )2) f (18) 

or from the temperature derivative of the average energy per edge (e) : 

c = d(e)/dT (19) 

In theory these formulas should give the same value, and the same should hold for measure- 
ments in the simulation, if the sampling is representative of the equilibrium. 

In Fig. [| the symbols represent the values of c obtained via formula ([181) . The errors 
are estimated by the blocking method. The continuous lines, on the other hand, come from 
fitting polynomials to the e-T plots in Fig. |] and then taking the deriviative with respect 
to T. We see that in the temperature range of the fits there is relatively good agreement 
between the two methods, thus suggesting that the simulation is equilibrated. 

The most notable feature of the plots are the peaks which reside near the temperature 
where m reaches 1. For N e = 45, 90 and 180 the peaks lie at T = 1.875, 1.375 and 1.125 
respectively 7 . For easier comparison, the curves for different N e are plotted together in Fig. 
|6]. The dotted line was added for clarity and connects the points by straight lines. 

This plot should be compared with Fig. |7|, which was obtained by taking the discrete 
forward derivative of the points in the m-T plot. The peaks of are slightly shifted 

relative to those of c. This may be an artifact resulting from the use of the forward derivative. 



7 It is not an accident that these values are close to 1. When choosing the couplings in the Hamiltonian, I 
tuned them to the onset of surface formation, while T was set to 1. 
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Figure 7: Forward derivative of measured values of m with respect to temperature T. 
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Figure 8: Susceptibility x versus temperature T. 



Figure [8] shows the results for the susceptibility x defined in eq. (||). x exhibits a sharp 
increase at the transition temperature. Unlike for a ferro magnet, there is no peak, however, 
since x stays more or less constant, as T is further increased. 

The salient properties of our data can be summarized as follows. Below a transition tem- 
perature T trans the surface observable m is constant and equals 1. At the temperatureT trans , 
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Figure 9: The fraction of edges with valence V e versus the temperature T for N e = 90. 

m starts to decrease. The rate of this decrease grows as N e increases. That is, the kink in 
the m-T curve becomes more pronounced as the system becomes larger. At the transition 
temperature T trans , we also see a peak in the specific heat c whose amplitude increases with 
system size. The peak in c is accompanied by a notable increase in the time At eq needed 
to equilibrate the system. A further indicator is the "susceptibility" x which shows a sharp 
rise at T = T trans . 

One can infer from this that for T < T trans the system is characterized by surface-like 
graphs for which m — 1. For T > T trans , these surfaces are replaced by more general graphs 
with m < 1. The observable m serves as a parameter to distinguish these two regimes. 

The T-dependence of the specific heat is similar to what is found in more conventional 
systems with phase transitions. It should be noted, however, that the peak of the specific 
heat shifts to lower temperatures, as the system size increases. Thus, it is possible that 
the transition temperature T trans goes to zero in the infinite size limit. Given that I have 
measured three sizes up to N e = 180 so far, it is premature to conclude whether this is the 
case or not. That is, one cannot tell at this stage if the system exhibits a phase transition at 
zero T or at finite T. Further remarks on this will be made in the conclusion of the paper. 



C. Surface properties 

This subsection deals in more detail with the surface states found at low temperature. 
What kind of surfaces appear? What defects do they have? And how do they change when 
T is increased? 

To answer these questions I conducted two additional sets of simulations. The first kind 
of simulation was aimed at showing more precisely how the valence of edges evolves as the 
temperature is raised above the transition temperature. This led to the plots in Fig. [| For 
T < T trans almost all edges have valence V e = 2. At T = T trans the fraction of these edges 
starts to decrease until it drops to around 66 percent at T = 2.0. In parallel the fraction 
of boundary edges (V e = 1) rises from percent to about 28 percent at T = 2.0. The 
percentage of edges outside of triangles (V e = 0) increases as well and reaches 5 percent. It 
is clear from this that at T = 2.0 the graphs must have lost most of its resemblance with a 
surface, since a considerable fraction of edges form boundaries of triangles. 
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0.40 ± 0.02 


N cc 


1.2 ± 0.1 


2.6 ± 0.3 


v v 


6.7 ± 0.1 


5.3 ± 0.1 


R 


-0.8 ± 0.1 


3.0 ± 0.3 


N lcJ N lc 


0.00 ± 0.00 


0.36 ± 0.01 


orient ability 


0.00 ± 0.00 


0.07 ± 0.01 


XEulcr 


-6.1 ± 0.4 


-4.8 ± 0.4 



Table II: Surface properties for edge number N e = 
The table consists of three sections which refer to 
graph) and the largest component in the surface. 



N v # of vertices 

N e # of edges 

Ny e= o # of edges with valence V e = 

-^con # of conical singularities 

N cc # of connected components 
of surface 

V v (average) valence of vertices 

R (average) Ricci scalar 

N l e c # of edges in largest component 

^le H 1 °f boundary edges in 
largest component 

180 and temperatures T = 0.625 and T = 2.0. 
the graph, the surface (defined by the reduced 



The second type of simulation concerned defects, topology and curvature. The measure- 
ments proceeded in three steps. First, the conical singularities and the edges with zero 
valence are determined and counted. Then, these defects are removed, so that one obtains a 
reduced graph that is equivalent to a triangulation of a surface. This reduction to a surface 
is needed in order to have a well-defined notion of topology. In the third step, we measure 
various quantities related to topology, boundary size and curvature of the surface. 

The curvature is obtained by interpreting the triangulation as a piecewise-linear manifold 
whose edges have length 1. The deficit angle at a vertex is given by 

5 V = 2?r - V v n/3 (20) 

and the Ricci scalar is denned by 



Here, A v = V v A/3 is the area associated to the vertex, and A = a/3/4 is the area of the 
triangles of the triangulation (see e.g. f26||). 

Since these measurements involve a change in the graph (and saving graphs was not 
implemented in the code), I ran the simulation differently from the previous ones. Instead 
of a sequence start — > equilibrate —> measure —> wait —> measure . . . , I used a schedule that 
restarts the simulation after each measurement: 

start — > anneal — > equilibrate — > measure 

— > restart — > anneal — > equilibrate — > measure . . . 

A similar simulation schedule will also play a role in the next section, where I discuss 
equilibration at low temperature. 
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Figure 10: Acceptance rate R a versus step number t for T = 0.5 and T = 2.0. 

The results are given in table |T|. I ran the simulation for T = 0.625 and T = 2.0 to 
compare the graph properties below and above the critical temperature. At low temperature, 
the graph contains no edges with valence V e = 0. Eight percent of the vertices come with 
conical singularities. The typical surface consists of a single connected component with 
moderate curvature. The largest component of the surface has no boundary, it is non- 
orientable and on average given by a connected sum of 8 real projective spaces. (Recall that 
a non-orientable surface with Euler characteristic XEuicr = 2 — k is a connected sum of k real 
projective planes.) 

At the temerature T = 2.0 above the transition point, the number of edges outside 
triangles is 9 times larger than in the surface regime, and the fraction of conical singularities 
reaches 40 percent. While we can still define a reduced graph and hence a surface, it is no 
longer a surface in the usual sense of the word. In the largest component, 30 to 40 percent 
of the edges are boundary edges. 

Let me add that I also tried to measure the spectral dimension of the graphs. This attempt 
failed, however. When measuring the return probability of diffusion, there is typically a 
window in the diffusion time that can be used to determine the spectral dimension. For 
diffusion times outside this interval the measurement is hampered by discretization and 
finite volume effects. In the present case, the system size appears to be so small (N e < 180) 
that discrete and finite volume effects overlap, and no clear signal of the spectral dimension 
is visible. 



V. SUBTLETIES AT LOW TEMPERATURE 

The simulation of systems with complex energy structure and many degeneracies can 
pose considerable challenges, especially at low temperatures. Well-known examples for this 
are systems with quenched disorder such as spin-glasses and random field systems (see 



e.g. fl20H22]). The simulation may get trapped in certain configurations, and never find 
the bottleneck to exit and enter other regions. It may sample some of the ground states, 
but leave out others, and thus create biased results that are not representative of the true 
equilibrium. 
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Figure 11: Comparison of the standard and alternate simulation schedule: on the left, specific heat 
c versus temperature T for standard schedule (grey) and modified schedule (black); in the middle, 
plot of e vs. simulation time t; on the right, parameters of Monte Carlo simulation. 



The ground states of the present model form a very special subset of the entire configu- 
ration space. Therefore, it is not far-fetched to imagine that similar difficulties could occur 
here. A first indication of this problem can be seen in the very small acceptance rates at 
low temperature. While the acceptance is already low, in general, it drops to values below 
10~ 3 as we reach T = 0.5 (see Fig. [H|). Real-time visualization with the Ubigraph soft- 
ware revealed that at this temperature the accepted changes are predominantly 2-2 Pachner 
moves — moves that are not able to change the topology of a triangulation. 

The question is therefore the following. Can the simulation sample the low-energy states 
in a statistically correct way, or does it get stuck in a minimum energy state and stay there 
for the rest of the simulation? To obtain more information about this, I repeated some of the 
previous simulations with a different schedule. In this new schedule, a phase of annealing and 
equilibration is followed by a sequence of measurements, then the simulation restarts and 
the whole process is repeated. This procedure may yield a better sampling of ground states 
and shorten the time during which the algorithm is trapped. The precise parameters of the 
original and modified schedule are indicated in Fig. [TT]. During annealing the temperature 
starts at T = 4.0 and drops gradually to the target temperature. 

The two types of simulations agree reasonably well for T > 1.0. For the lowest tempera- 
tures, however, the alternate schedule leads to a much higher value of the specific heat. The 
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Figure 12: Autocorrelation function a(t) of e for different temperatures and N e = 180. 



reason for this difference becomes clear when we inspect the measured data for the energy 
per edge (see Fig. [II]). In the standard case, the fluctuations exhibit a fixed pattern for the 
length of the entire simulation. In the alternate case, the same kind of pattern appears, 
however, only for the duration of a single measurement period. At the beginning of the fol- 
lowing measurement sequence, the pattern jumps to a new energy and remains there until 
the next restart. 

This suggests the following interpretation. The standard schedule gets stuck in one set 
of states and stays there for the rest of the simulation. The alternate schedule gets trapped 
as well, but only for a short period and each time for a different energy. This sampling of 
different states results in a too large variance and hence in a wrong value for the specific 
heat. Within this sample, however, there are also states that have lower energy than the 
states found by the conventional scheme. Thus, there is something wrong and right in both 
methods. The original simulation yields the correct specific heat, but the measured energy 
may be too high, if the algorithm freezes in a set of states above the true ground state. 

The last point to be addressed is the autocorrelation of our data (see Fig. |12|). For 
temperatures greater than T trans , the autocorrelation drops exponentially to zero and I de- 
termined corresponding correlation times At cor for the error analysis. For T < T trans , on the 
other hand, the data are, in many cases, autocorrelated on time scales of the order of t MC , 
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so I would not be able to correct this other than by going to much larger simulation times. 
The plots show furthermore that my estimates for the equilibration time were too optimistic 
at low temperatures. 

VI. CONCLUSION 

Let me summarize the results. I have defined a graph model that produces "space from no 
space". Two-dimensional triangulations (space) arise as ground states from a much larger 
set of configurations (no space). To measure the extent to which a graph is a surface I 
introduced an observable called m. Monte Carlo simulations show that below a transition 
temperature T trans the observable m is close or equal to 1, and the graphs form triangulations 
of surfaces 8 . For temperatures T > T trans , one has m < 1 and the surfaces are replaced by 
more general graphs. More detailed measurements reveal that the surfaces are typically 
connected and non-orientable. As we go above the transition temperature, the graphs lose 
the surface structure, since the number of defects increases and many of the edges become 
boundaries of triangles. 

The temperature dependence of m and of the specific heat c suggest that there is a 
phase transition between the regime with and without surfaces. However, the transition 
temperature T trans drifts, as the system size is increased, so it is possible that T trans goes 
to zero in the infinite-size limit. In this case, there would be no phase transition at finite 
temperature, and one would instead speak of a zero temperature phase transition. Since the 
system was only simulated for very few and small sizes, it is too early to tell which scenario 
is realized. 

What are future directions for this work? Obviously, one should try to extend the model 
to larger sizes and higher, more realistic dimensions. The key technical question in this 
will be whether one is able to make the algorithm efficient enough to sustain equilibrium at 
low temperature. It may be necessary to invent more sophisticated simulation schedules in 
order to deal with larger edge numbers and higher-dimensional simplices. Another remedy 
could be to use a different Hamiltonian that is more "flexible" in its choice of ground states. 
By this I mean that ground states would not only be given by graphs that are literally 2d 
triangulations, but also by more general graphs that are effectively 2-dimensional (e.g. in 
the spectral sense). Perhaps one could define a "spectral Hamiltonian" that achieves this. 

Furthermore, one should determine conclusively whether there is a phase transition, and 
if it is at finite or zero temperature (in the infinite-size limit). A phase transition at finite 
temperature is desirable in the sense that it would make the behaviour near the transition 
universal and reduce the amount of fine-tuning needed to achieve this state. In the case of 
a zero temperature phase transition, the dimensionless couplings would have to be sent to 
infinity as the system size goes to infinity in order to maintain a transition at finite T. 

Since the paper focused on the technical aspects, I would like to conclude with a number 
of conceptual remarks. An obvious question pertains to the role of time. If this is a statistical 
model of space, there must be a time with respect to which the system equilibrates. If that 
is the case, it would seem that, in the model's world, time is of a different nature than space, 
since space can dissolve, while there is always time for the system to thermalize. 

Let me give a possible interpretation of the model that clarifies this point. Imagine that 



modulo conical singularities 
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we had a consistent quantum theory of gravity that describes spacetime as being built from 
elementary cells or building blocks. Assume also that from a canonical point of view this 
spacetime corresponds to blocks or "atoms" of space that evolve in a discrete, fluctuating 
time. If the theory is a viable candidate theory of quantum gravity, there should be a 
regime where classical continuum spacetime is recovered from the collective properties of 
these atoms. Like in conventional quantum theories, one might furthermore expect that 
there is a statistical approximation of the theory, where the 4d path integral is replaced by 
a statistical sum over 3d configurations. Our model could be regarded as such a statistical 
approximation to a quantum theory of spacetime (in this case 3d spacetime). The edges of 
the graphs are the atoms of space and their interactions are encoded by the Hamiltonian 9 . 

This line of thought suggests also a connection with group field theory and loop quantum 
gravity. The quanta of these theories are triangles and tetrahedra and interact to form 
pseudo-manifolds p7|-^9||. If we were able to develop a statistical theory of group field 
theory and loop quantum gravity, the resulting models may have a similar structure as the 
model defined here 10 . 

Apart from this, our system could be also useful as a concrete and simple toy model 
that helps to imagine what a physics beyond spacetime could be; a physics, that is, in 
which spacetime is no longer absolute and instead only a special state among other states 
of the system. This dethroning of spacetime is, in my view, a logical continuation of the 
development of theoretical physics. It is also a radical step that challenges accustomed ways 
of thinking about physics. Let me highlight this by describing three thought experiments 
that could be simulated with our model. 

There is a Id version of the Hamiltonian that results in the formation of Id chains at low 
temperature. Imagine that we combine the Id Hamiltonian with the 2d Hamiltonian used 
in this paper. Then, one could vary the relative strength of the Id and 2d couplings and 
and observe transitions between a 1-dimensional and 2-dimensional phase of the system. 

As a second example, consider the surface regime of the model and suppose that we 
encode a picture on the surface by affixing additional labels to edges. Next let us raise 
the temperature until the surface disappears completely. Then, we go back to the original 
temperature and a surface will form again. The picture, however, will be scrambled, since 
the arrangement of edges is completely reshuffled. One can view this as an instance of 
spontaneous symmetry breaking. 

Another version of this argument explains nicely the intimate connection between space 
and Hamiltonian (or spacetime and Lagrangian). Suppose that we extend the system by 
adding spins to the triangles and let these spins interact by the usual nearest-neighbour 
Hamiltonian. Assume furthermore that the energy scale of the spin field is small enough to 
not affect the dynamics of the graphs. At low enough temperatures, the graphs freeze in 
a surface state and the system becomes effectively an Ising model on a fixed lattice. Note, 
however, that it is the graph that determines which spin neighbours which other spin and 
hence which spins couple to each other. In this sense, the graph (or space) is part of the 
spin system's coupling constants. By the same token, if space is a variable, the coupling 
constants are variables as well. 



9 Clearly, this is just a heuristic consideration and leaves out many questions. What effect would fluctuations 

of time have? How can Lorentz invariance be maintained? 
See |30| for another possible link between loop quantum gravity and statistical physics. 
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